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(57) A fluid property monitor includes a transducer 
assembly (T, 26') to impart multiple frequency energy to 
a conduit in one or more modes and to receive (T, RE- 
CEIVER) resonant frequency energy from the conduit. 
The resonant frequency energy is responsive to the im- 
parted energy, the conduit and a fluid in the conduit. The 
monitor can also include: a frequency signal generator 
connected to cause multiple frequency energy to be 



transferred to a conduit having a fluid to be monitored; 
and a spectral analysis signal processor connected to 
receive and process electrical signals generated in re- 
sponse to vibrations propagated through the conduit 
and the fluid in the conduit in response to transferred 
multiple frequency energy. Particular implementations 
can be adapted as a densitometer, a coherent flow de- 
tector, and other particular fluid parameter detectors. 
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Description 

[0001] This invention relates generally to fluid proper- 
ty monitors and monitoring methods. In particular, the 
invention provides a fluid property monitor and method 
that use multiple frequency inputs in one or more vibra- 
tional modes to produce multiple resonant frequency 
outputs from which one or more fluid characteristics are 
determined. A particular application is to determine den- 
sity of a cement slurry for use in an oil or gas well, but 
other applications are also contemplated. 
[0002] Cement slurry is mixed on rigs prior to delivery 
into wells. A cement slurry typically needs to have a par- 
ticular density or density range to control downhole 
pressures during cementing. Slurry density is now typi- 
cally monitored with nuclear densitometers; however, 
nuclear sources are forbidden in some countries. Also, 
density readings are sometimes erratic for foam ce- 
menLs at the wellhead after adding nitrogen. Other den- 
sitometers with fewer regulatory requirements (e.g. Co- 
riolis devices) may disrupt flow and suffer from erosion 
by fluids. Still others (e.g. high frequency acoustic de- 
vices) may yield only small depths of investigation into 
the cross section of a fluid, which is unacceptable for 
fluids having radial gradients in pipes, for example. In 
other applications, density, compressibility and viscosity 
measurements are needed to monitor various flowing 
fluids, such as stimulation fluids and produced fluids 
downhole in wells. 

[0003] A fluid monitor and monitoring method mitigat- 
ing or overcoming the aforementioned shortcomings are 
needed. 

[0004] In one aspect, the invention provides a fluid 
property monitor, comprising a transducer assembly to 
impart multiple frequency energy to a conduit in multiple 
modes, and to receive resonant frequency energy from 
the conduit, wherein the resonant frequency energy is 
responsive to the imparted energy, the conduit and a 
fluid in the conduit. 

[0005] The monitor of the invention preferably imparts 
vibrations to a distinct system of a conduit and fluid, from 
which the monitor detects responses resulting from in- 
duced deformations of the conduit and the effects on the 
fluid in the conduit. 

[0006] The vibrations used in the present invention re- 
late to the modes of vibration of cylindrical shells. Al- 
though standing wave patterns for various such modes 
may be well known, the use and implementation of mul- 
tiple modes and/or multiple frequency response in de- 
termining fluid characteristics as called for by the 
present invention are new. For a flexural mode, fluids in 
the bore of shells (to a first approximation) merely add 
mass loading, resulting in a decrease of resonant fre- 
quency. Radial mode resonance is affected by speed of 
sound which depends on density and compressibility. 
Torsional mode response couples primarily through the 
shear viscosity. The Q, the ratio of the center frequency 
divided by the difference between the upper and lower 



frequencies having amplitudes 3 dB smaller than the 
peak amplitude, of any of the resonances relates to en- 
ergy dissipation, which relates to fluid shear viscosity, 
fluid bulk viscosity, other fluid dissipation mechanisms 

5 and structural energy loss mechanisms. Q's of the mul- 
timodes can be solved by multiple regression to esti- 
mate the contributions of each energy loss mechanism. 
[0007] Such a fluid monitor and method may be used 
to measure various parameters, including density, com- 

10 pressibiiity and viscosity (and preferably measured from 
the same sample volume). Accuracy and response 
times of the invention preferably equal or exceed those 
of prior devices for monitoring corresponding fluids and 
parameters. The presence of multiple resonances per- 

15 mits the use of several center frequencies to improve 
the statistical accuracy of the determined characteristic. 
The invention has no nuclear sources, and it avoids de- 
vice erosion. 

[0008] The present invention applies to mixing and 
20 delivery of cement slurry, for example. As another ex- 
ample, the invention may also apply to evaluating stim- 
ulation fluids during delivery for enhanced reservoir pro- 
duction. Additionally, the invention maybe used for com- 
bining with fluid velocity measurements to develop mass 
25 flow meters for various applications, including perma- 
nent downhole sensors. 

[0009] Particular embodiments are described in fol- 
lowing sections, but some other advantages of these in- 
clude the following. 

30 [0010] Fixed positions between transducers providing 
vibrational energy give relatively stable resonant fre- 
quencies even if the locations of end clamps supporting 
the test region conduit change slightly. The use of force 
couples avoids lateral forces on the end clamps that oc- 

35 cur with unbalanced translational forces (forces on end 
clamps may cause errors in estimating true resonant fre- 
quencies, which frequencies are used in determining 
the desired fluid characteristic(s)). 
[0011] Separate transducers for transmitting and re- 

40 ceiving provide greater noise rejection than combined 
transmitter/receiver approaches; however, combined 
transducer embodiments are encompassed in broader 
aspects of the present invention. 
[0012] Using a relatively large number of transmitters 

45 permits large amplitude vibrations to be generated. The 
dependence of amplitude on frequency gives the non- 
linear behavior of fluid modulii which relates to the 
amount of gas in liquid. The invention, therefore, may 
be useful for estimating the amount of gas in liquids. 

so [0013] One embodiment permits transducers to be 
retrofitted on pipes that are already installed. The inven- 
tion, therefore, can be designed with no interruption of 
the flow contours. Interruption of flow contours can 
cause erosion or improper mixing of multi-phase fluids. 

55 [0014] In measuring at least density and using at least 
the flexural mode, the invention monitors the entire liq- 
uid cross section and therefore gives the average den- 
sity of the total flow. 
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[001 5] Another advantage is non-intrusive monitoring 
of density or other parameters. The radial flow profile in 
the test region conduit is unaffected by the device. This 
device requires no curvature of the conduit axis as in 
curved Coriolis flow meters. The maximum conduit di- 
ameter is only limited by the length of the monitor appa- 
ratus and the flexural wavelength. The larger the diam- 
eter, the longer the monitor apparatus. Active length of 
the activated conduit can be reduced relative to flexural 
mode length by using only one or more of the radial 
modes in implementing the present invention. 
[0016] One definition of the present invention is as a 
fluid property monitor comprising a transducer assem- 
bly to impart multiple frequency energy to a conduit in 
multiple modes (e.g., flexural modes and/or radial 
modes), and to receive resonant frequency energy from 
the conduit, wherein the resonant frequency energy is 
responsive to the imparted energy, the conduit and a 
fluid in the conduiL. In a particular embodiment, a radial 
mode includes a radial hoop mode and a radial oval 
mode. 

[001 7] The fluid property monitor can also be defined 
as comprising: a frequency signal generator connected 
to cause multiple frequency energy to be transferred to 
a conduit having a fluid to be monitored; and a spectral 
analysis signal processor connected to receive and 
process electrical signals generated in response to vi- 
brations propagated through the conduit and the fluid in 
the conduit in response to transferred multiple frequen- 
cy energy. 

[0018] In a particular implementation the present in- 
vention provides a densitometer that comprises: at least 
four transducers disposed circumferentially relative to a 
location along a length of a conduit connected in-line in 
a flow path for a fluid to be monitored; a multiple fre- 
quency signal, generator to provide one or more drive 
signals forthe transducers to impart multiple frequency 
energy to the conduit; and a controller to connect the 
one or more drive signals to the transducers such that 
the transducers are operated to drive the conduit in any 
of a plurality of modes including a multiple frequency 
flexural mode and a multiple frequency radial mode. 
[0019] In another implementation the present inven- 
tion provides a coherent flow detector comprising first 
and second fluid property monitors disposed at respec- 
tive first and second locations along a conduit for a fluid 
to be monitored. The first fluid property moniLor includes 
a first transducer assembly to impart first multiple fre- 
quency energy to the conduit at the first location in mul- 
tiple modes, including at least a flexural mode and a ra- 
dial mode ; and to receive first resonant frequency ener- 
gy from the conduit, wherein the first resonant frequency 
energy is responsive to the imparted first multiple fre- 
quency energy, the conduit and a fluid in the conduit. 
The second fluid property monitor includes a second 
transducer assembly to impart second multiple frequen- 
cy energy to the conduit at the second location in mul- 
tiple modes, including at least a flexural mode and a ra- 



dial mode, and to receive second resonant frequency 
energy from the conduit, wherein the second resonant 
frequency energy is responsive to the imparted second 
multiple frequency energy, the conduit and the fluid in 

5 the conduit. 

[0020] A method of monitoring a fluid in accordance 
with the present invention comprises: imparting multiple 
frequency vibration-inducing energy to a conduit and flu- 
id system; and sensing a plurality of frequency signals 

10 from the conduit and fluid system responsive to at least 
part of the imparted multiple frequency vibration-induc- 
ing energy. This method can also comprise: determining 
at least one characteristic of the conduit and fluid system 
in response to the sensed plurality of frequency signals. 

15 in a particular implementation, determining at least one 
characteristic includes using an averaging calculation 
with frequencies identified from the sensed plurality of 
frequency signals. 

[0021] A method of monitoring a fluid in accordance 
20 with the present invention can also be defined as com- 
prising: imparting multiple frequency energy at a first lo- 
cation of a conduit for a fluid to be monitored; imparting 
multiple frequency energy at a second location of the 
conduit; and detecting vibrational motion at a third loca- 
ls tion of the conduit, wherein the third location is between 
the first and second locations. 

[0022] Another definition of the method of monitoring 
a fluid comprises: flowing a fluid through a conduit; driv- 
ing at least one of one or more transducers disposed 

30 adjacent the conduit to impart energy to the conduit to 
deform the conduit in a flexural mode: driving at least 
one of the one or more transducers to impart energy to 
the conduit to deform the conduit in a radial hoop mode; 
driving at least one of the one or more transducers to 

35 impart energy to the conduit to deform the conduit in a 
radial oval mode; and generating signals, with at least 
one of the one or more transducers, in response to the 
driving steps. 

[0023] Therefore, from the foregoing, it is a general 
40 object of the present invention to provide a novel and 
improved vibrational forced mode fluid property monitor 
and method. Other and further objects, features and ad- 
vantages of the present invention will be readily appar- 
ent to those skilled in the art when the following descrip- 
45 tion of the preferred embodiments is read in conjunction 
with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

50 [0024] FIG. 1 is a schematic illustration of a conduit 
with which a plurality of transducers are associated and 
driven to impart flexural mode vibrational energy to the 
conduit and a fluid inside the conduit. 
[0025] FIG. 2 is a schematic illustration of the conduit 

55 with which the plurality of transducers are associated 
but now driven to impart radial mode vibrational energy 
to the conduit and a fluid inside the conduit. 
[0026] FIG. 3 is a schematic cross-sectional illustra- 
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tion showing a circumferential group of four transducers 
driven to provide a radial hoop mode. 
[0027] FIG. 4 is a schematic cross-sectional illustra- 
tion showing the circumferential group of four transduc- 
ers driven to provide a radial oval mode. 5 
[0028] FIGS. 5A-5F represent sequential stages with- 
in a cycle of the radial oval mode, wherein FIG. 5A rep- 
resents a starting time, FIG. 5B depicts a maximum ver- 
tical elongation of the conduit, FIG. 5C depicts reversed 
drive polarity, FIG. 5D represents the conduit being driv- 10 
en back through its neutral shape, FIG. 5E depicts a 
maximum horizontal elongation of the conduit, and FIG. 
5F illustrates the next reversed drive polarity to start the 
next cycle upon returning to the neutral shape of FIG. 
5A. 15 
[0029] FIG. 6 is a schematic illustration of a conduit 
with which three sets of four circumferentially spaced 
transducers each are associated such that the transduc- 
ers can be driven to impart flexural, radial hoop, and ra- 
dial oval modes of vibrational energy to the conduit and 20 
a fluid inside the conduit. 

[0030] FIG. 7 is a schematic illustration of a moment 
arm assembly connected to a conduit to impart torsional 
mode vibrational energy to the conduit. 
[0031] FIG. 8 is a schematic illustration of moment -25 
arm assemblies connected to a conduitto.impart flexural 
mode vibrational energy to the conduit. 
[0032] FIG. 9 is a schematic representation of a single 
pair of transducers driven to impart a flexural mode to a 
conduit. 30 
[0033] FIG. 10 is a schematic representation of the 
single pair of transducers driven to impart a radial mode 
to the conduit. 

[0034] FIG. 1 1 is a schematic representation of a sin- 
gle combination transmitter/receiver transducer used to 35 
impart a flexural mode and a radial mode to the conduit 
and to receive responsive signals. 
[0035] FIG. 12 is a schematic and block representa- 
tion of a coherent flow detector of the present invention 
comprising two fluid property monitors of the present in- 40 
vention. 

[0036] FIG. 13 illustrates a particular implementation 
having five pairs of transducers mounted on the conduit 
and driven to impart flexural mode vibrations. 
[0037] FIG. 14 illustrates the particular implementa- 45 
tion driven to impart radial mode vibrations. 
[0038] FIG. 15 is a block diagram of transmitting and 
receiving portions of the fluid property monitor of the 
present invention to be described as particularly corre- 
lated to the particular implementation of FIGS. 13 and 50 
14. 

[0039] FIG. 16 is an elevational view of a particular 
implementation of a transducer for the embodiment of 
FIGS. 13 and 14. 

[0040] FIG. 17 is a block diagram of transmitting cir : 55 
cuitry to drive the piezoelectric transducers of the imple- 
mentation of FIGS. 13 and 14. 

[0041 ] FIG. 1 8 is a schematic diagram of a mode con- 
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troller circuit connected between the transmitting circuit- 
ry of FIG. 1 7 and the drive transducers. 
[0042] FIG. 1 9 is a schematic circuit diagram of a sig- 
nal conditioning circuit connected to a receiver trans- 
ducer of the implementation of FIGS. 1 3 and 14. 
[0043] FIG. 20 shows the results from averaging cal- 
culated densities for each oval mode resonant frequen- 
cy in tests using the identified fluids in a static state. 
[0044] FIG. 21 is a schematic diagram of a test flow 
system incorporating the implementation of FIGS. 13 
and 14. 

[0045] FIG. 22 shows graphs of resonant frequencies 
detected during a test in which flow rate of the fluid 
changed in the system of FIG. 21 . 
[0046] FIG. 23 shows graphs of resonant signal am- 
plitudes across the flow rate range. 
[0047] FIG . 24 shows a frequency domain data acqui- 
sition example for normalized densitometer response 
over a broad frequency range used to describe a math- 
ematical analytical technique in the signal processing of 
one implementation of the present invention. 
[0048] FIG . 25 shows a definition for a nine point mov- 
ing average used in the mathematical analytical tech- 
nique to smooth received resonant signals. 
[0049] FIG. 26 provides a densitometer frequency re- 
sponse plot with multi-layer Fibonacci searches em- 
ployed to find center resonant frequencies in the math- 
ematical analytical technique. 

[0050] FIG. 27 illustrates part of a log normalized res- 
onant signal for which half-power amplitude bandwidth 
is to be found in the mathematical analytical technique. 

DETAILED DESCRIPTION OF THE INVENTION 

[0051] The vibrational forced mode fluid property 
monitor and method of the present invention is for de- 
termining one or more various parameters for various 
fluids. Examples of such characteristics include density, 
compressibility, viscosity, flow velocity (including slip ve- 
locity, such as with two monitors used as a coherent flow 
detector), and foaming; and examples of fluids include 
cement slurries (including foamed slurries), stimulation 
fluids (e.g., fracturing fluids, acids), and drilling muds. 
The fluids can be static or flowing, and the fluids can be 
Newtonian or non-Newtonian and homogeneous or 
non-homogeneous. Use can be at various locations 
along a flow stream (e.g., before or after nitrogen injec- 
tion for a foamed cement slurry) and with respect to var- 
ious locations within the fluid (e.g., at the center, at the 
outer boundary, or an average across a cross section). 
To simplify the explanation of the invention, however, 
reference in the following description will primarily be 
made to monitoring density of a cement slurry; however, 
such specifics are not to be taken as limiting broader 
aspects of the invention. 

[0052] The modes to which the present invention per- 
tains include flexural, radial, and torsional, with the first 
two of preferred interest. At present, the most preferred 
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embodiment uses a combination of fiexural mode, radial 
hoop mode, and radial oval mode; however, any of the 
modes can be used alone or in combination with one or 
more of the other modes. When the monitored fluid has 
multiple phases (e.g., gas and liquid), fiexural, hoop, 5 
and oval modes have different responses dependent 
upon the spatial distribution of the fluid phases. The fiex- 
ural mode is more sensitive to fluid phases at the center 
of the stream than are the two mentioned radial modes. 
The hoop mode changes the volume of the flow conduit 10 
more than does the oval mode; therefore, the hoop 
mode is more sensitive to compressibility effects due to 
the presence of gases in liquids than are the fiexural or 
oval modes. 

[0053] Within the modes, multiple frequencies are is 
used. Multiple frequency mechanical energy is applied 
or imparted to the conduit containing a fluid under test 
such that multiple frequency vibrations are induced in 
the conduit and fluid. Frequencies throughout or across 
the range of imparted frequencies are monitored and the 20 
relatively larger amplitude frequencies (typically the res- 
onant frequencies) are used to determine the one or 
more parameters or characteristics (e.g., density). That 
is. the present invention senses for responsive signals 
related to at least some of the imparted frequencies, 25 
namely, for the resonant frequencies produced by cor- 
responding input frequencies in the preferred embodi- 
ments. 

[0054] For resonant frequencies associated with fiex- 
ural modes, important factors include pipe structure, 30 
transducer mass, mass inside the pipe, and mass out- 
side the pipe. For unchanging interior fluids, fiexural 
mode frequencies can be calibrated to measure density 
(for example) of fluids outside the pipe. Previous fiexural 
mode approaches position drivers for modes governed 35 
by the pipe length. Resonances based on pipe length 
have errors due to greater dependence on rigidity of end 
clamps. The forced mode approach of preferred embod- 
iments of the present invention reduces end effects by 
having a central receiver with multiple drivers between 40 
the receiver and either pipe end. Driver axial locations 
are selected to enhance signals at receivers and mini- 
mize effects of end clamps. In FIG. 1 transducer pairs 
2, 4 are equidistantly spaced from transducer pair 6. To 
provide the fiexural mode, the pairs are driven as indi- 45 
cated by the respective arrows and then reversed in a 
cyclical manner. Although imparted from respective lo- 
cations along the conduit, the cumulative effect applies 
vibratory energy to the conduit 8 as indicated by axial 
and boundary lines 10, 12, 14 (the illustrated displace- 50 
ments are exaggerated relative to the actual sizes and 
distances). 

[0055] Resonant frequencies of radial modes depend 
on pipe structure, transducer mass, fluid inside and fluid 
outside thepipe, andfluid speed of sound. FIG. 2 shows 55 
the transducer pairs from FIG. 1 used to create a radial 
mode. The radial mode created by driving the transduc- 
ers 2, 4, 6 as shown by the arrows in FIG. 2 and then 



reversing the drive in a cyclical manner propagates by 
deflecting the wall as shown in FIG. 2 by boundary lines 
1 2 : 14. The wall deflection at receiver transducer pair 6 
is proportional to the deflection at drivertransducerpairs 
2, 4, but delayed in time. Driving two sets of drivers at 
equal symmetric distances relative to the radial plane at 
the receivers gives increased signal amplitude. Care 
must be taken if more sets of drivers at other axial dis- 
tances from the receivers are used. The different prop- 
agation distances for other axial distances can give de- 
structive rather than constructive interference at the re- 
ceiver. The phase mismatch depends on speeds of 
sound in the fluid and in the conduit wall. 
[0056] Two distinct radial modes are preferably gen- 
erated with different electrical phases for the transduc- 
ers. Considerfourtransducers T equally spaced circum- 
ferentially around the pipe as shown in FIGS. 3 and 4. 
Driving all four transducers T in phase causes uniform 
expansion and contraction of the circumference- to de- 
fine a radial hoop mode (FIG. 3). Driving one opposing 
transducer pair in phase and the other opposing trans- 
ducer pair 180 degrees out of phase results in a radial 
oval mode (FIG. 4, and see FIGS. 5A-5F for an illustra- 
tion of a cycle of the oval mode deflection of the conduit). 
These modes have fundamentals and overtones/that 
give two distinct sets of resonant frequencies. Process- 
ing the sets independently gives greater accuracy of the 
density (or other characteristic) measurement. 
[0057] FIG. 6 represents an embodiment with two 
sets 1 6, 1 8 of four circumf erentially spaced driver trans- 
ducers each. These are equally spaced axially from a 
centrally located set of four circumferentially spaced re- 
ceiver transducers 20. By adjusting the relative electri- 
cal phasing of the individual driver and receiver trans- 
ducers, one can vibrate the pipe 22 in fiexural, hoop or 
oval mode. Preferably each of these modes is used and 
responsive frequencies detected to determine the de- 
sired parameter (e.g., density). Should end effects arise 
despite the symmetrical arrangement depicted, they 
would be more troublesome for the fiexural mode than 
for the radial modes; however, the combination of flex- 
ural and radial modes permits using the radial modes to 
compensate the fiexural mode for end effects. 
[0058] When the transducers of FIG. 6 are driven as 
in FIG. 3 to create a hoop mode, the central axis of the 
conduit 22 is not moved as in the fiexural mode, making 
the system relatively insensitive to density of material at 
the center. This mode has greater sensitivity to density 
of material at the wall. Because the area of the circle 
expands and contracts, the volume of the pipe changes 
and this mode is sensitive to changes in the net com- 
pressibility of non-homogeneous material anywhere in 
the pipe. When the transducers of FIG. 6 are driven as 
in FIGS. 4 and 5 to create an oval mode, the central axis 
of the pipe 22 is not moved. The area of the oval chang- 
es little ; reducing the effects of compressibility. 
[0059] The foregoing are preferred modes of the 
present invention; however, another mode included in 
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the present invention with regard to modes used in com- 
bination with each other is the torsion mode. For reso- 
nant frequencies associated with torsion modes., impor- 
tant factors include pipe structure, transducer mass, 
mass and shear viscosity of fluid inside the pipe, mass 
and shear viscosity of fluid outside the pipe, fluid com- 
pressibility, and transducer location. FIG. 7 shows a mo- 
ment arm arrangement, including transducers T and 
clamped moment arm members 26, that imparts torsion- 
al motion to pipe 24. One of these arrangements placed 
at the axial center of the pipe will create torsional dis- 
placement of the pipe wall . The receiver must sense tor- 
sional displacements. Transverse oriented accelerom- 
eters, transverse poled piezoelectrics or torsion arm 
transducers are acceptable. Torsion modes couple 
weakly through viscosity and so a non-torsional mode 
may obscure effects of the weak coupling; therefore, a 
torsion driver should be designed to avoid non-torsional 
excitation modes and a receiver should be designed to 
reject small non-torsion modes created by torsion driv- 
ers. 

[0060] The moment arm arrangement can also be 
used to induce the flexural mode as illustrated in FIG. 
8. The forced mode generated in this device can be de- 
veloped by applying a set of equal and opposite couples 
at equal distances from the center of the device. The 
moments can be created by mounting the transducers 
on moment arms 26'attached to the pipe as shown in 
FIG. 8. The transducers oscillate in the axial direction of 
the pipe and thus do not co ntribute any force in the trans- 
verse direction relative to the pipe. The moment impart- 
ed by the oscillation of the transducers and the flanges 
bends the pipe while allowing vibration in preferred nat- 
ural modes. Such a vibration is represented in FIG. 8 by 
center line 28. The pipe is forced to observe a change 
in bending moment at the clamping flanges: The pipe 
response to the driving frequency is a strong function of 
the pipe contents, since the pipe vibrates in its natural 
modes. To obtain an antinode at the receivers, the driv- 
ers must produce pipe wall motions that are symmetric 
relative to a radial plane at the centered receivers. The 
two moments shown in FIG. 8 acting in opposite direc- 
tions and located on opposite sides of the center remove 
any asymmetric mode shapes and allow only symmetric 
modes to exist. The centered antinode of symmetric 
modes can be easily detected by the receivers mounted 
at the center of the pipe. The centered nodes of asym- 
metric modes, in contrast, have no transverse deflec- 
tions at the center of the pipe, and hence give no signals 
in the receivers. A receiver similar to that for the moment 
arm drivers in FIG. 7 is used to detect centered nodes 
of asymmetric modes. 

[0061] FIG. 8 shows transducers on moment arms. 
An alternative design includes a flange that allows res- 
onance at desired frequencies using either normal or bi- 
morph piezoelectric elements. 

[0062] The foregoing has been described with refer- 
ence to more than three transducers. In this particular 



context, the fluid property monitor of the present inven- 
tion comprises a receiver disposed relative to a length 
of conduit to detect vibrational energy of a location of 
the conduit (e.g., a central location as illustrated in 
5 FIGS. 1, 2, 6, and 8). This fluid property monitor also 
comprises a first transmitter to impart vibrational energy 
to the conduit. The first transmitter is disposed relative 
to the conduit at a location spaced along the length of 
the conduit in a first direction from the receiver. This fluid 
to property monitor still further comprises a second trans- 
mitter to impart vibrational energy to the conduit. The 
second transmitter is disposed relative to the conduit at 
a location spaced along the length of the conduit in a 
second direction from the receiver such that the receiver 
f5 is between the first and second transmitters. For exam- 
ples of these transmitters, see FIGS. 1, 2, 6, and 8. 
[0063] Although the foregoing represents presently 
preferred embodiments of the present invention, one 
can create and receive aL leasl flexural and radial modes 
20 with a single transducer pair 30a, 30b mounted relative 
to conduit 32 as illustrated in FIGS. 9 and 10; however, 
. the interpretation of the respective modes may be diffi- 
cult if the range of density (or other measured parame- 
ter) causes overlapping of the response ranges for the 
25 respective modes; Furthermore, even a single combi- 
nation transmitter/receiver transducer 34 can be used 
as illustrated adjacent conduit 36 in FIG. 11 . This capa- 
bility of generating many modes with a single transducer 
results because many modes cause radial movement 
30 of the pipe surface, and a transducer receiving signals 
related to the radial movement detects the modes. 
Transducers generating forces in the radial direction ex- 
cite the modes if driven at the resonant frequency (trans- 
ducers placed at nodes or other locations giving no ra- 
35 dial wall motion cannot receive or generate the corre- 
sponding modes). In view of these options, what is 
needed at a minimum for the fluid property monitor of 
the present invention is a transducer assembly to impart 
multiple frequency energy to a conduit in one or more 
40 modes (e.g., flexural and/or radial modes), and to re- 
ceive resonant frequency energy from the conduit, 
wherein the resonant frequency energy is responsive to 
the imparted energy, the conduit and a fluid in the con- 
duit. In terms of a method, this includes imparting mul- 
45 tiple frequency vibration-inducing energy to a conduit 
and fluid system, and sensing a plurality of frequency 
signals from the conduit and fluid system responsive to 
at least part of the imparted multiple frequency vibration- 
inducing energy. In response to the sensed plurality of 
50 frequency signals, at least one characteristic of the con- 
duit and fluid system is determined. As will be explained 
below, determining at least one characteristic includes 
using an averaging calculation with frequencies identi- 
fied from the sensed plurality of frequency signals. In 
55 the embodiments of FIGS. 1,2,6 and 8, for example, 
multiple frequency energy is imparted at respective lo- 
cations along the conduit and resultant vibrational mo- 
tion is detected at another location of the conduit in be- 
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tween the first-mentioned locations. 
[0064] The present invention also provides two or 
more fluid property monitors as described above to cre- 
ate a coherent flow detector that can indicate the speed 
at which inhomogeneities move relative to the monitors 5 
for multiphase fluids. Such a coherent flow detector 
comprises: a first fluid property monitor disposed at a 
first location along a conduit for a fluid to be monitored, 
and a second fluid property monitor disposed at a sec- 
ond location along the conduit. In one implementation, 10 
the first fluid property monitor includes a first transducer 
assembly to impart first multiple frequency energy to the 
conduit at the first location in multiple modes, including 
at least a flexural mode and a radia! mode. The first 
transducer assembly is also to receive first resonant fre- *5 
quency energy from the conduit, wherein the first reso- 
nant frequency energy is responsive to the imparted first 
multiple frequency energy, the conduit and a fluid in the 
conduit. In such implementation the second fluid prop- 
erty monitor includes a second transducer assembly to 20 
impart second multiple frequency energy to the conduit 
at the second location in multiple modes, including at 
least a flexural mode and a radial mode. The second 
transducer assembly is also to receive second resonant 
frequency energy from the conduit, wherein the second 25 
resonant frequency energy is responsive to the impart- 
ed second multiple frequency energy, the conduit and 
the fluid in the conduit. This is schematically illustrated 
in FIG. 12. 

30 

A Particular Implementation 

[0065] A specific fluid property monitor which has 
been built and used in internal experimentation referred 
to hereinbelow is represented in FIGS. 13 and 14. 35 
Transducers 40a, 40b, 42a, 42b, 44a : 44b, 46a, 46b are 
placed along conduit or pipe 48 to give balanced forces 
at least in flexural mode. For flexural mode, two pairs of 
outer transducers (40, 46) produce forces in the oppo- 
site direction to two pairs of inner drivers (42, 44) as rep- 40 
resented in FIG. 13. In FIG, 13, the four pairs 40, 42, 44, 
46 produce two couples, one couple on either side of 
centered receiver transducers 50a, 50b. Each couple 
has no net translational forces, only rotational forces. 
This feature means that end clamps experience super- 4 $ 
positions of waves from two couples. Spacing between 
couples and between each couple and pipe ends are 
selected to minimize end effects. 

[0066] For radia! mode, at least one of the pairs 40, 
42, 44, 46 is used, wherein the respective transducers 50 
of the pair(s) are cyclically driven towards and away 
from each other as depicted in FIG. 14 for pair 42. If only 
one pair is used, the other driving transducers are qui- 
escent. The behavior in FIG. 1 4 can be achieved by an 
electrical system that allows the electrical phase of each 55 
transmitter and each receiver to be actuated either in 
phase or 1 80° out of phase with respect to an oscillator 
driven at the selected modal resonance frequency. 



[0067] When the conduit and transducer assembly of 
FIGS. 13 and 14 is to be used, the driver transducers 
are connected to receive one or more multiple frequen- 
cy, multiple phase drive signals provided from a multiple 
frequency signal generator through a mode controller 
as indicated in FIG. 15. The controller connects one or 
more drive signals to one or more of the drive transduc- 
ers so thatthetransducers impart the multiple frequency 
energy to the conduit in the appropriate mode(s). The 
electrical outputs from the receiver transducers re- 
sponding to what has happened to the imparted energy 
in the conduit and fluid system are provided to a frequen- 
cy analyzer from which one or more measured param- 
eters is determined as also indicated in FIG. 15. The 
frequency analyzer is a signal processor that provides 
spectral analysis of the received signals, particularly 
resonant frequency energy detection responsive to the 
imparted energy, the conduit and the fluid in the conduit. 
The features indicated in FIG. 1 5, including the conduit 
48 (FIGS. 13 and 14) and the transducers 40, 42, 44, 
46. 50, will next be described in more detail for this par- 
ticular implementation of a fluid property monitor of the 
present invention. 

[0068] The conduit 48 can be a section of pipe already 
in place and containing a flowing fluid. Transducers;ean 
be clamped to or otherwise operatively associated with 
the outside of the pipe. The system can be calibrated by 
flowing materials with known density through the pipe. 
Alternatively, a flanged section can be installed which 
has been p re-calibrated. The conduit 48 can be a stand- 
alone test chamber or pipe segment or other suitable 
fluid container. 

[0069] Because large amplitudes are desired to inter- 
pret the nonlinear response of the liquid modulii.the pre- 
ferred embodiment has distributed transducers to mini- 
mize pipe distortion. The large forces for a single driver 
could cause local distortion of the pipe and so this is not 
preferred, but, as stated above, the present invention 
can be implemented with a single transducer. 
[0070] In the implementation of FIGS. 13 and 14, to 
avoid coupling changes with aging, the preferred ap- 
proach is to weld nuts to fasten transducers to the pipe; 
however, other attachment techniques can be used (e. 
g., clamping), and non-contacting transducers can be 
used. An example of a non-contacting transducer is a 
permanent magnet attached to the pipe and a drive coil 
mounted on a stationary frame. The coil is energized al 
the driving frequency. 

[0071] Referring to FIG. 16, each of the transducers 
used in the particular construction related to FIGS. 13 
and 14 mentioned above is a studded transducer that 
attaches to a respective nut 60 welded to the conduit 
48. A cylindrical threaded post 62 extending from both 
sides of a stainless steel base 64 has a locking nut 65, 
a tensioning nut 66, and a Bellville spring washer 68 
mounted on it to press a piezoelectric member 70 to- 
wards the base 64. The washer 68 acts against a tung- 
sten washer 72 which provides inertia! backing for the 
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piezoelectric 70. The washer 72 lies against an insulat- 
ing washer 74 disposed adjacent a negative electrode 
76 positioned against one side of the piezoelectric mem- 
ber 70. A positive electrode 78 and an insulating washer 
80 are between the other side of the piezoelectric mem- 5 
ber 70 and the base 64. The portion of the post 62 ex- 
tending" on this side of the base 64 passes axially 
through the elements as shown in FIG. 1 6. An insulating 
sleeve 82 is between this portion of the post 62 and the 
elements up to the tensioning nut 66. The other portion 10 
of the post 62 extends from the opposite face of the base 
64 to screw into the nut 60. 

[0072] The insulatingwashers 74, 80 on the transduc- 
ers functioning as receivers give common mode rejec- 
tion of electrical crosstalk from the drivers. If added re- 15 
jection is needed, a stack of piezoelectrics can be used 
with the top side of the top element having the same 
electrical polarity as the bottom side of the bottom ele- 
ment. The greater the number of elements in the stack, 
the better the rejection of capacitive electrical crosstalk. 20 
[0073] If an embodiment using this particular con- 
struction is retrofitted on an existing pipe that precludes 
welded nuts, the stainless steel base 64 can be made 
without the lower (as viewed in FIG. 16) stud and con- 
toured to the outside diameter of the pipe. The tension- 25 
ing and locking nuts 65, 66 can be replaced with a clamp 
that encircles the pipe. The clamp pulls together two 
transducers mounted on diametrically opposite sides of 
the pipe. The transducer placement would be the same 
as that in FIGS. 13 or 14, but the pairs of transducers 30 
would be secured by encircling clamps that pressed in- 
ward on th e Bellville washers instead of by studs in weld- 
ed nuts. 

[0074] If the positive (+) receiver in FIG. 13 has posi- 
tive polarity on the outside of the piezoelectric stack and 35 
the negative (-) receiver has negative polarity, then" 
these receivers attach in parallel to the electronic am- 
plifier in the frequency analyzer represented in FIG. 15. 
The parallel configuration further reduces capacitive 
crosstalk. Because receivers have stacks of positive 40 
and negative polarities, radial mode vibrations are re- 
jected. Sensitivity to flexural mode vibrations, however, 
is unaffected. This configuration reduces capacitive 
crosstalk and mechanical detection of undesired radial 
modes. (If switching from flexural to radial mode is de- 45 
sired, receivertransducers are individually connected to 
separate amplifiers.) 

[0075] If translational mechanical vibrations interfere 
with the flexural motion, a second pair of receivers can 
be mounted at an odd number of wavelengths from the so 
illustrated receiver pair 50. The locations of the positive 
and negative receivers for this second pair are inter- 
changed diametrically. All receivers are attached elec- 
trically in parallel. The four receivers should have elec- 
trical rejection of capacitive crosstalk and mechanical 55 
rejection of radial modes and translational modes. This 
wiring configuration allows cancellation without the 
noise of separate amplifiers. Amplifier noise comes after 



the electrical subtraction. It also avoids the dynamic 
range problems amplifiers have in subtracting large sig- 
nals that are almost equal. This four receiver configura- 
tion may be needed for environments with large me- 
chanical shocks but may be unnecessarily expensive for 
normal conditions. If the amplifier dynamic ranges are 
inadequate for subtraction, separate sets of receivers 
for flexural, for hoop and for oval mode can be used to 
obtain subtraction before amplification. 

Transmitting Circuitry 

[0076] For the transducers used as drivers in the par- 
ticular implementation being described with reference 
to FIGS. 13-15, the transducer drive circuitry provides 
the multiple frequency, multiple phase energy signals 
and mode control indicated in FIG. 1 5. Circuitry for gen- 
erating the drive signals will next be described with ref- 
erence to FIG. 17, and then an implementation for'the 
mode controller will be described with reference to FIG. 
18. In this implementation power for the circuitry is uni- 
polar (e.g., +12VDC from a vehicle battery), but split 
power (+/-15VDC in a specific implementation) is used 
for piezoelectric transducer drive. 

[0077] FIG. 1 7 depicts a particular implementation for 
generating sine wave signals used to drive the trans- 
ducers associated with the conduit in which the fluid to 
be tested is contained or through which the fluid to be 
tested flows. A triangle wave is provided from a triangle 
wave generator 90. An alternative control signal can be 
provided from a programmed microprocessor circuit, or 
microcontroller, 92 (e.g., Z-Worid BL1700-based cir- 
cuit). Selection of the control signal is made by the mi- 
crocontroller 92 via its operation of a solenoid switch 94 
as apparent from FIG. 17. The selected control signal is 
input to a voltage-to-frequency converter 96, the fre- 
quency output from which is divided to a desired level 
by binary clock dividers 98. The output of the circuit 98 
drives a sine wave generator (with bandpass filter) 1 00. 
The output of the sine wave generator 100 is provided 
to the mode controller. This output from the sine wave 
generator 1 00 is also provided to the microcontroller 92 
for monitoring and synchronization. 
[0078] In the particular embodiment, the foregoing cir- 
cuitry is designed and operated to generate drive sig- 
nals that vibrate the conduit 48 within the range between 
about 100 hertz and about 20 kilohertz (kHz), and more 
preferably between about 3 kHz and about 14 kHz. Such 
frequencies can be obtained by providing one or more 
signals having a frequency that is sweptthrough the par- 
ticular range used or by providing multiple individual sig- 
nals each having a respective frequency within the 
range or by providing a composite multiple frequency 
signal (e.g., a multiple frequency pulse) or a combina- 
tion of the foregoing. 

[0079] Depicted in FIG. 1 8 is an embodiment of the 
mode controller for the particular implementation being 
described with reference to the general block diagram 
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of FIG. 15. The output signal from the sine wave gener- 
ator 1 00 is filtered to remove direct current (DC) biases 
and attenuate any high frequency noise above 14kHz 
for the aforementioned preferred range. In this context, 
filtering can be purely passive because gain control is 
not necessary. This output signal is split eight ways (i. 
e., provided in parallel to eight inputs) and routed 
through parallel inverting/non-inverting amplifier sets 
102a. 102b ...102h (e.g., model HA-2400 devices in a 
specific implementation) that offer gain and phase con- 
trol. Outputs from these amplifier sets are coupled to- 
gether and routed to power operational amplifier circuits 
104a : 104b, ...104h (e.g., model OPA544 devices in a 
specific implementation). Final drive signals are then 
routed through transformers 106a, 106b, ...106h to the 
eight piezoelectric driver transducers 40, 42, 44, 46 of 
FIGS. 13 and 14. Each transformer 106 preferably ad- 
equately matches both source (power op amp 1 04) and 
destination (transducer) impedances, offers minimal 
distortion over the specified frequency range, and han- 
dles the estimated six-watt (W) power load with reason- 
able size, weight, and cost. Phase, and thus mode, con- 
trol occurs by the microcontroller 92 controlling which 
amplifier within a set 1 02 is operational to process and 
output the input drive signal. Different combinations of 
in-phase and/or 180 out-of-phase signals can be pro- 
vided. 

[0080] Advantages of this approach include the fol-_ 
lowing. For phase control, switching between amplifiers 
1 02 occurs at low power ratherthan at high power which 
should reduce noise tendencies. Supply current de- 
mands are met on the low voltage side, which is more 
easily accomplished than on the high voltage side. 
[0081] Although the foregoing is preferred, an alter- 
native is to use high voltage operational amplifiers for 
final transducer drive, powered by a high voltage power 
supply. This would eliminate the audio transformers 1 06 
of FIG. 18, but then needed high voltage linear power 
supplies rated at roughly 0.1 amp (A) rms per transducer 
(0.8 A rms total) at 60-70 volts (V) rms would be very 
large, expensive items with a large power transformer 
included. 

Receiving Circuitry 

[0082] Turning now to the receiving side of the pres- 
ently described particular implementation for the em- 
bodiment of FIG. 15, a high impedance isolation ampli- 
fier can be used for piezoelectric receiver signal condi- 
tioning. That is, the output of a receiver transducer con- 
nects to a signal conditioning amplifier, and one type of 
such amplifier is a high impedance isolation amplifier. 
For example, input impedance on the order of 80 meg- 
aohms provides crystal isolation, with signal voltage am- 
plification available through a resistor pair ratio. Another 
type of amplifier, and one that is presently preferred, is 
a low input impedance current amplifier. An example of 
a receiver amplifier circuit is shown in FIG. 19. Other 



conditioning components include a full-wave rectifier, 
differential amplifier, and low-pass integrator. The signal 
output from the receiver circuit is processed for frequen- 
cy analysis which leads to measurement of the param- 

5 eter(s). Such analysis and interpretation into parameter 
or characteristic measurement is described further be- 
low with regard to processing techniques, which can be 
implemented using the microcontroller 92 of FIG. 17. 
The resultant measurement can be output by any suit- 

10 able means; non-limiting examples include a display 
screen (e.g., a touch screen which also allows user in- 
put), a 4-20 mA interface to provide data in this conven- 
tional format (as in a radioactive densitometer), and a 
serial interface (such as for digital data transfer). 

15 

Test System 

[0083] A test system for a fluid property monitor of the 
• type shown in FIGS. 3 to 6 was suspended vertically in 

20 the laboratory with a valve on the bottom. The tube was 
filled with the valve closed. Resonant frequencies were 
measured, and the tube was emptied by opening the 
valve. The electrical driving oscillator was controlled 
manually to find resonant frequencies for a variety of liq- 

25 uids. For the oval mode in FIG. 5, the results in FJG.- 20 
were obtained. 

[0084] A second fluid property monitor 108 (FIG. 21) 
was used to gauge resonant frequency and bandwidth 
in a flowing system. This monitor used theflexural mode 

30 shown in FIG. 1 3. FIG. 21 shows a centrifugal pump and 
an alternating current (AC) induction electric motor 110 
with adjustable pulleys allowing for three different pump 
speed ranges. Flow rate adjustment within each pump/ 
motor speed range was achieved by adjusting the loop 

35 restriction via a valve 112. This allowed the motor to run 
at constant speed. Connection hoses and a reservoir 
114 were also included as apparent from the drawing. 
To complete the flow loop, a Halliburton turbine flow me- 
ter 1 1 6 was acquired, and various pipe fittings were pur- 

40 chased and assembled. Water (specific gravity 1 .0) was 
pumped from the reservoir 1 1 4 by the motor/pump 110, 
through the fluid property monitor 108 and the flow me- 
ter 116, and back into the reservoir 114. 
[0085] Two dependent variables were of primary in- 

45 terest with flow rate the independent variable. A 0.5-in. 
Halliburton turbine flow meter was used with conversion 
constant 875 pulses per gallon. Flow rate (gallons per 
minute) was computed using this constant along with 
frequency measurements taken during flow tests. Res- 

50 onant frequencies near 6200 Hz, 9000 Hz, and 12000 
Hz were measured by centering their resonant spike on 
an oscilloscope graticule and measuring the associated 
voltage-controlled oscillator (VCO) signal voltage. Sys- 
tem calibration data was used to convert VCO ampli- 

55 tudes to frequency. Some measurement uncertainty 
was inherent, but of primary interest was to see if reso- 
nant frequencies had any flow rate dependence. Also 
measured were the actual resonant spike voltage am- 
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plitudes. Of interest was to see how amplitudes were 
maintained through a series of tests. 
[0086] Flow rates from 2 gpm (gallons per minute) to 
1 6 gpm were achieved. The graphs of FIG. 22 show that 
frequency for the most part held constant through the 5 
flow rate range. Aeration became more significant as 
flow rates increased, possibly due to leakage through a 
pump seal. One expects frequency to increase with per- 
cent aeration (overall fluid density decreases), but only 
the 1 2000 Hz trace exhibited that behavior. Average fre- 10 
quencies were 6228, 9076, and 12053 Hz. 
[0087] FIG. 23 shows much variation in signal ampli- 
tude through the flow range. This was first believed to 
be exclusively a product of aeration . Above 1 2 gpm , aer- 
ation was significant and all signal amplitudes decayed is 
greatly. At approximately 5.9 gpm, however, all signal 
amplitudes were attenuated through a narrow band. 
Possible causes for this behavior are at present un- 
known; however flow velociLy at 5.9 gpm was approxi- 
mately 6.8 in./sec, and it is possible that 6.8 in. /sec was 20 
coincident with some critical flow velocity. 
[0088] The foregoing illustrates that in the present in- 
vention the frequency response for a given spacing of 
transducers is determined by sweeping driving frequen- 
cy and observing center frequencies of resonances: A 25 
number of resonance peaks will be observed, repre- 
senting various reverberation modes of the entire struc- 
ture. Peaks with largest amplitudes and largest mechan- 
ical Q's are observed for various calibration fluids from 
air to the most dense fluid to be measured. The reso- 30 
nances having the most stable amplitudes and Q's for 
the range of calibration densities should be used for in- 
terpreting density (or other fluid characteristic). The per- 
centage change of center frequency versus character- 
istic is relatively stable and can be calibrated for each 35 
resonance peak. For improved accuracy, the character- 
istic's measurement can be obtained from a statistical 
fit of several resonant peaks. A statistical estimate 
based on the center frequencies of standing waves may 
give greater accuracy than resonance peaks individual- 40 
Iv- 

[0089] Increasing vibrational amplitude eventually 
creates nonlinear response which gives changing cent- 
er frequency versus amplitude. By measuring change 
of frequency versus amplitude, the nonlinear character- 45 
istics of the pipe and drivers can be compensated. Ad- 
ditional nonlinear behavior of the system due to adding 
liquid can then be determined relative to the response 
of calibration liquids such as water. The calibration liq- 
uids should be single phase with density comparable to so 
the liquids to be characterized. This calibration proce- 
dure is useful for characterizing two phase liquids which 
have a nonlinear modulus and may be useful for identi- 
fying small percentages of gas in liquid. 
[0090] Monitoring multiple resonant modes as de- 55 
scribed above can help reduce the effects of resonanc- 
es from mechanical noises. Mechanical noises will be 
out of synchronization with the sweep rate of the driver 



transducers and therefore have incoherent effects on 
the multiple resonant frequencies. A regression fit to 
multiple center frequencies statistically reduces effects 
of mechanical noise at resonance frequencies. 

Processing Techniques 

[0091 ] The foregoing gives a sample of the use of the 
present invention; however, an actual implementation 
takes the respective detected resonant frequencies, 
amplitudes and associated bandwidths (preferably the 
3dB bandwidth, or Q) and converts them to the desired 
parameters (e.g., density or any of the other character- 
istics previously mentioned). Two techniques for doing 
this will next be described, but it is contemplated that 
other techniques (including a combination of empirical 
and mathematical) can be used and remain within the 
scope of the present invention. The two techniques to 
be described here are the use of an empirical neural nel- " 
work and the use of explicit mathematical algorithms. 

Empirical Technique 

[0092] An equation for estimating density, for exam- 
ple, from the i th resonance frequency is: 



Densityj = Aj / (Frequency^ 4- B } 

[0093] The preferred determination of the density is 
obtained by averaging the results for all relevant reso- 
nate frequencies for each respective mode (i.e., aver- 
age of the above A/(Frequencyj) 2 +B| for all the identified 
resonant frequencies for the respective mode), and av- 
eraging between or among two or more modes can be 
used if they relate to the same characteristics to be de- 
termined. For example, a flexural mode test can be per- 
formed to determine its correlation to a particular char- 
acteristic, density, for example. This is done with one or 
more known fluids, and from this the A and B values can 
be determined. The density for an unknown fluid can 
then be determined by using the A and B values with 
the resonant frequencies identified in a flexural mode 
test of the unknown fluid. An average is calculated using 
the individual values derived for each resonant frequen- 
cy. The average can also include values derived from 
any other relevant mode; therefore, the present inven- 
tion includes averaging within each mode (which alone 
may be used) and/or averaging using two or more 
modes. When multiple modes are used, the initial cal- 
culation for each mode is performed in the same manner 
as for the flexural mode mentioned above. FIG. 20 
shows the results for averaging the calculated densities 
for each oval mode resonant. frequency in a test using 
the identified fluids in a static state as described above. 
The vertical pipe test system was used on the fluids list- 
ed on the graph. The straight line manually fitted to the 
points is within ± 0.015 gm/cc (gram per cubic centim- 
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eter) of the measured data points. 
[0094] The values for the Aj's and Bj's of the above 
equation can be obtained by using calibration fluids with 
known densities and measured resonant frequencies 
for each mode of interest (flexural, hoop, oval torsional). 5 
These values of calibration constants are then used to 
calculate the relevant characteristics of unknown fluids. 
The Aj's and B/s are different functions of compressibil- 
ity and viscosity. They are determined using fluids hav- 
ing known densities, compressibilities, and viscosities. 10 
Center frequencies and Q's responsive to the multiple 
frequency vibrations applied to the fluid test system are 
determined. The frequencies and Q's are dependent da- 
ta input for a neural network computer. The independent 
data are the fluid densities, compressabilities. and vis- is 
cosities forthe calibration fluids. During a calibration run 
the neural network establishes predictive relations. 
From the predictive relations, the programmed compu- 
ter then calculates densily and other parameters for un- 
known fluids to which multiple frequency vibrations are 20 
applied in the test system. 

Mathematical Technique 

[0095] The following describes a mathematical ap- 25 
proach in correlating detected response to actual pa- 
rameter. Initially it is noted that frequency response 
(magnitude) curves are multi-modal about each reso- 
nance, that is, several relative maxima appear. Center 
frequency selection (absolute maxima) is complicated 30 
by presence of relative maxima. Described below is a 
technique to locate absolute maxima for any resonant 
frequency band, leading to proper center frequency 
identification. This methodology is then adapted to 
bandwidth limit identification. 35 

Center Frequency Identification 

[0096] FIG. 24 provides a frequency domain data ac- 
quisition example that shows normalized densitometer 40 
response over a broad frequency range. Integrated re- 
ceiver response can be somewhat noisy with multi-mo- 
dal peaks. The challenge is to scan data, properly pick 
the absolute peak, then identify the corresponding res- 
onant frequency. With multi-mode data, misidentifica- 45 
tion of the absolute maximum peak must be avoided. 
[0097] Some smoothing is initially used to condition 
all piezoelectric receiver data. In one test, a nine point 
moving average (defined in FIG. 25) was used. All nine 
weights (wtj) add to unity (1 .0). Multi-mode traces still 50 
existed after smoothing, but some point-to-point noisi- 
ness was reduced. A Fibonacci search method was then 
used to locate peak amplitudes. 

[0098] FIG. 26 provides a densitometer frequency re- 
sponse plot with multi-layer Fibonacci searches em- 55 
ployed. All three peak amplitudes were properly identi- 
fied, leading to proper center frequency identification. 
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Bandwidth identification 

[0099] Many of the problems associated with center 
frequency identification also apply to bandwidth identi- 
fication. Detection of upper and lower bandwidth fre- 
quency limits fy and f L is useful for interpreting energy 
dissipation of the fluid, which relates to viscosity. 
[0100] Referring to FIG. 27, assume that center fre- 
quency f c = f[ifc] with index 'ifc' having been identified. 
The half-power amplitude lies at approximately -0.3 on 
a log-normalized scale. Lower bandwidth limit index 'ifP 
can be readily identified, from which f L = f[ifl] is achieved. 
Note, however, that several crossings of the -0.3 thresh- 
old can occur when noisy data exists. In this case, the 
best solution is to identify all -0.3 crossings so that an 
overall, 'average' bandwidth limit frequency can be de- 
termined. 

[0101] A modified Fibonacci-type search can be ap- 
plied. The center frequency index ifc can be used as the 
reference point. Lower bandwidth limit data index 'if I! can 
be identified by searching data to the left (lower indices) 
and upper bandwidth limit data index 'ifu' can be identi- 
fied by searching data to the right (upper indices). When 
searching each peak for bandwidth limits, the iogrnor- 
malized amplitude must be considered. The data set 
must be either re-normalized or bandwidth limits must 
be adjusted. 

Conversion of Frequency to Density 

[0102] Frequency uncertainty Af n arises from two 
sources. One uncertainty source is born from signal dis- 
cretization that is unavoidable. The other source relates 
to center frequency selection, but this source of uncer- 
tainty can be avoided or reduced by using the mathe- 
matical technique referred to above. Once output fre- 
quencies have been determined, however, they can be 
converted to measurements of the respective charac- 
teristics using the formula and neural network approach . 
referred to above. In general, however, each identified 
resonant frequency is correlated to the respective char- 
acteristic and an average of all the values is calculated 
to define the characteristic measurement that is to be 
output or otherwise used. This is performed within each 
excitation mode and across all respective modes used 
as desired. 

[0103] Accordingly, one aspect of the present inven- 
tion provides a method of monitoring a fluid comprising: 
imparting mechanical waves into a fluid to be monitored, 
wherein the mechanical waves have frequencies within 
a predetermined range; generating alternating current 
electric signals in response to transmission character- 
istics of the fluid relative to the imparted mechanical 
waves; and indicating a characteristic of the fluid in re- 
sponse at least to identified frequencies within the alter- 
nating current electric signals. In the preferred embodi- 
ments described above, the identified frequencies are 
resonant frequencies of a system which includes theflu- 
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id. A particular implementation of the method more spe- 
cifically comprises: flowing a fluid through a conduit; 
driving at least one of one or more transducers disposed 
adjacent the conduit to impart energy to the conduit to 
deform the conduit in a flexural mode; driving at least 5 
one of the one or more transducers to impart energy to 
the conduit to deform the conduit in a radial hoop mode; 
driving at least one of the one or more transducers to 
impart energy to the conduit to deform the conduit in a 
radial oval mode; and generating signals, with at least 10. 
one of the one or more transducers, in response to the 
driving steps. This can further comprise indicating at 
least one characteristic of the flowing fluid in response 
to the generated signals. 

[0104] Thus, the present invention is well adapted to 15 
carry out the objects and attain the ends and advantag- 
es mentioned above as well as those inherent therein. 
While preferred embodiments of the invention have 
been described for the purpose of this disclosure, 
changes in the construction and arrangement of parts 20 
and the performance of steps can be made by those 
skilled in the art, which changes are encompassed with- 
in the spirit of this invention as defined by the appended 
claims. 
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Claims 

1. A fluid property monitor comprising a transducer 
assembly to impart multiple frequency energy to a 30 
conduit in multiple modes, and to receive resonant 
frequency energy from the conduit, wherein the res- 
onant frequency energy is responsive to the impart- 
ed energy, the conduit and a fluid in the conduit. 

35 

2. A monitor according to claim 1 , further comprising: 
a signal generator connected to the transducer as- 
sembly to cause at least part of the transducer as- 
sembly to vibrate at multiple frequencies: and a sig- 
nal processor connected to the transducer assem- *o 
bly to receive from at least part of the transducer 
assembly electrical signals responsive to the reso- 
nant frequency energy. 



or more drive signals for the at least four transduc- 
ers; and a controller connected to the at least four 
transducers and the signal generator to connect the 
one or more drive signals to one or more of the at 
least four transducers so that the transducers im- 
part the multiple frequency energy to the conduit in 
the multiple modes. 

7. A fluid property monitor, which comprises a fre- 
quency signal generator connected to cause multi- 
ple frequency energy to be transferred to a conduit 
having a fluid to be monitored; and a spectral anal- 
ysis signal processor connected to receive and 
process electrical signals generated in response to 
vibrations propagated through the conduit and the 
fluid in the conduit in response to transferred multi- 
ple frequency energy. 

8. A monitor according to claim 7. wherein the multiple 
frequency energy is acoustic energy within the" 
range of 1 00 hertz to 20,000 hertz. 

9. A monitor according to claim 7 or 8, wherein the 
transducer means includes a flexural mode and a 
radial mode, in response to the frequency signal 
generator for generating and providing the electrical 
signals to the spectral analysis signal processor. 

1 0. A fluid property monitor which comprises a receiver 
disposed relative to a length of a conduit to detect 
vibrational motion of a location of the conduit, 
wherein the conduit is for a fluid to be monitored; a 
first transmitter to impart vibrational energy, prefer- 
ably of multiple frequencies, to the conduit, the first 
transmitter disposed relative to the conduit at a lo- 
cation spaced along the length of the conduit in a 
first direction from the receiver; and a second trans- 
mitter to impart vibrational energy to the conduit, the 
second transmitter disposed relative to the conduit 
at a location spaced along the length of the conduit 
in a second direction from the receiver such that the 
receiver is between the first and second transmit- 
ters. 



3. A monitor according to claim 1 or 2, wherein the ra- 
dial mode includes a hoop mode and an oval mode. 

4. A monitor according to claim 1 , 2 or 3, wherein the 
transducer assembly includes a single combination 
transmitter/receiver transducer, or a plurality of 
transducers disposed along a length of the conduit. 

5. A monitor according to claim 4, wherein the trans- 
ducer assembly includes at least four transducers 
disposed circumferentially relative to the conduit. 

6. A monitor as defined in claim 5 : further comprising 
a multiple frequency signal generator to provide one 



45 11. A monitor according to claim 1 0, further comprising 
means for operating the first and second transmit- 
ters such that the conduit is deformed in multiple 
modes, including a flexural mode and a radial 
mode, in response to the imparted vibrational ener- 

50 gy. 

12. A fluid property monitor, which comprises at least 
three transducer groups disposed at respective lon- 
gitudinal locations along a conduit connected in-line 
55 . j n a flow path for a fluid to be monitored; a multiple 
frequency signal generator; a controller connecting 
the signal generator to the transducer groups such 
that a selectable output from the signal generator is 
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provided to one or more transducers of a selected 
one or more of the transducer groups; and a signal 
processor connected to at least one transducer of 
at least one of the transducer groups to receive and 
process electrical signals generated by the at least 5 
one transducer to which the signal processor is con- 
nected in response to a selected output from the 
signal generator driving the one or more transduc- 
ers of the selected one or more of the transducer 
groups. 10 

13. A densitometer which comprises at least four trans- 
ducers disposed circumferentially relative to a loca- 
tion along a length of a conduit connected in-line in 

a flow path for a fluid to be monitored; a multiple '5 
frequency signal generator to provide one or more 
drive signals for the transducers to impart multiple 
frequency energy to the conduit; and a controller to 
connect the one or more drive-signals to the trans- 
ducers such that the transducers are operated to 20 
drive the conduit in any of a plurality of modes in- 
cluding a multiple frequency flexura! mode and a 
multiple frequency radial mode, said multiple fre- 
quency radial mode preferably including a multiple 
frequency hoop mode and a multiple frequency oval 25 
mode. 

14. A densitometer according to claim 1 3, further com- 
prising a signal processor connected to at least one 

of the transducers to receive and process electrical 30 
signals generated by the at least one transducer, to 
which the signal processor is connected, in re- 
sponse to vibrational energy in the conduit resulting 
from the multiple frequency mode in which the con- 
duit is driven. 35 

1 5. A coherent flow detector which comprises a first flu- 
id property monitor disposed at a first location along 
a conduit, the first fluid property monitor including a 
first transducer assembly to impart first multiple f re- 40 
quency energy to the conduit at the first location in 
multiple modes, including at least a flexural mode 
and a radial mode, and to receive first resonant fre- 
quency energy from the conduit, wherein the first 
resonant frequency energy is responsive to the im- ^5 
parted first multiple frequency energy, the conduit 
and a fluid in a conduit; and a second fluid property 
monitor disposed at a second location along the 
conduit, the second fluid property monitor including 

a second transducer assembly to impart second so 
multiple frequency energy to the conduit at the sec- 
ond location in multiple modes, including at least a 
flexural mode and a radial mode, and to receive 
second resonant frequency energy from the con- 
duit wherein the second resonant frequency ener- 55 
gy is responsive to the imparted second multiple fre- 
quency energy, the conduit and the fluid in the con- 
duit. 



16. A method of monitoring a fluid, which method com- 
prises imparting multiple frequency vibration-induc- 
ing energy to a conduit and fluid system; preferably 
by driving the conduit and fluid system in at least 
one of a flexural mode and a radial mode; and sens- 
ing a plurality of frequency signals from the conduit 
and fluid system responsive to at least part of the 
imparted multiple frequency vibration-inducing en- 
ergy. 

17. A method according to claim 1 6, further comprising 
determining at least one characteristic of the con- 
duit and fluid system in response to the sensed plu- 
rality of frequency signals, preferably by using an 
averaging calculation with frequencies identified 
from the sensed plurality of frequency signals. 

18. A method as defined in claim 16 or 1 7, wherein driv- 
ing the conduit and fluid system in -a- radial mode 
includes driving the conduit and fluid system in at 
least one of a radial hoop mode and a radial oval 
mode. 

19. A method of monitoring a fluid which method com- 
prises imparting multiple frequency energy at.a first 
location of a conduit for a fluid to be monitored: im- 
parting multiple frequency energy at a second loca- 
tion of the conduit; and detecting vibrational motion 
at a third location of the conduit, wherein the third 
location is between the first and second locations. 

20. A method according to claim 1 9, further comprising 
controlling the imparting of multiple frequency en- 
ergy at the first and second locations so that the 
conduit is deformed in multiple modes, including a 
flexural mode and a radial mode, in response to the 
imparted energy, said method preferably further 
comprising determining at least one characteristic 
(e.g. density) of the fluid in response to the detected 
vibrational motion. 

21 . A method of monitoring a fluid, which method com- 
prises driving transmittertransducers located along 
a conduit symmetrically relative to a receiver trans- 
ducer connected to the conduit: and determining at 
least one characteristic, preferably the density, of a 
fluid in the conduit in response to a plurality of fre- 
quency signals from the receiver transducer re- 
sponsive to driving the transmittertransducers. 

22. A method of monitoring a fluid, which method com- 
prises imparting mechanical waves into a fluid to be 
monitored, wherein the mechanical waves have fre- 
quencies within a predetermined range; generating 
alternating current electric signals in response to 
transmission characteristics of the fluid relative to 
the imparted mechanical waves; and indicating a 
characteristic, preferably the density, of the fluid in 
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response at least to identified frequencies within the 
alternating current electric signals, said identified 
frequencies preferably being resonant frequencies 
of a system which includes the fluid. 

5 

23. A method of monitoring a fluid, which method com- 
prises flowing a fluid through a conduit; driving at 
least one of one or more transducers disposed ad- 
jacent the conduit to impart energy to the conduit to 
deform the conduit in a flexural mode: driving at 10 
least one of the one or more transducers to impart 
energy to the conduit to deform the conduit in a ra- 
dial hoop mode; driving at least one of the one or 
more transducers to impart energy to the conduit to 
deform the conduit in a radial oval mode; and gen- is 
erating signals, with at least one of the one or more 
transducers, in response to the driving steps; and 
preferably indicating at least one characteristic, 
preferably the density, of Lhe flowing fluid in re- 
sponse to the generated signals, preferably by iden- 20 
tifying resonant frequencies in the generated sig- 
nals for each mode; and determining the at least 
one characteristic in response to an averaging cal- 
culation using the identified resonant frequencies. 
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